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In this work we have performed for the first time a systematic analysis of the disso- 
ciation and fission pathways of neutral, singly and doubly charged alanine dipeptide 
ions with the aim to identify the fission mechanism and the most probable frag- 
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t^ . mentation channels of these type of molecules. We demonstrate the importance of 

rearrangement of the molecule structure during the fission process. This rearrange- 
ment may include transition to another isomer or a quasi-molecular state before 
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f->i ■ actual separation of the daughter fragments begins. 
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Amino acids are building blocks for proteins. Recently, it became possible to study 
experimentally fragments of proteins, i.e. chains of amino acids, in a gas phase with the 



use of the MALDI mass spectroscopy 
acid chains began also only recently 




Abj/nitio theoretical investigations of amino 
10| and are still in its infancy. In 
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this work we have performed for the first time a systematic analysis of the dissociation and 
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fission pathways of neutral, singly and doubly charged alanine dipeptide ions with the aim to 
identify the fission mechanism and the most probable fragmentation channels of these type 
of molecules. We demonstrate the importance of rearrangement of the molecule structure 
during the fission process. This rearrangement may include transition to another isomer or 
a quasi-molecular state before actual separation of the daughter fragments begins. 

II. THEORETICAL METHODS 

Our exploration of the energetics and the pathways of the fragmentation process is based 
on the density-functional theory (DFT) and the post-Hartree-Fock many-body perturbation 
theory accounting for all electrons in the system. Within the DFT one has to solve the 
Kohn-Sham equations, which read as (see e.g. |lOl.lll|): 



y + Uions + Vh + Vxc 1 V'i = ^iV'i, (1) 

where the first term represents the kinetic energy of the z-th electron, and f/jons describes its 
attraction to the ions in the cluster, Vu is the Hartree part of the interelectronic interaction: 

VHir-)= fj^dr', (2) 

/ ly ly I 

and p{r') is the electron density: 

N 
u=l 

where Vxc is the local exchange-correlation potential, ipi are the electronic orbitals and A^ is 
the number of electrons in the cluster. 

The exchange-correlation potential is defined as the functional derivative of the exchange- 
correlation energy functional: 

_ 5E,c[p] ,. 

^^' ~ bp{f) ' ^^^ 

The approximate functionals employed by DFT methods partition the exchange- 
correlation energy into two parts, referred to as exchange and correlation parts. Both parts 
are the functionals of the electron density, which can be of two distinct types: either lo- 
cal functional depending on only the electron density p or gradient-corrected functionals 
depending on both p and its gradient, Vp. In literature, there is a variety of exchange 



correlation functionals. In our work we use the Becke's three parameter gradient-corrected 
exchange functional with the gradient-corrected correlation functional of Lee, Yang and Parr 
(B3LYP) |l2, ll^, llJ| • We utilize the standard 6-31G(d) basis set to expand the electronic 
orbitals -04. 

To explore the energetic and the pathways of the fragmentation process we introduced the 
procedure which implies the calculation of the multidimensional potential energy surface. 
We increase the distance between the two fragments and calculate the multidimensional 
potential energy surface on each step. Then we identify local minima on this surface for 
a given distance between the fragments, and so find the fragmentation pathways of the 
molecule. For our calculations we utilized the Gaussian 03 software package |l5Lllfil|. 

III. RESULTS OF CALCULATION 

In our work we have investigated the fission process of the two types of the alanine 
dipeptide. The only difference between them is in the N— ending of the molecule, which 
could be terminated by the —NH2 or —NH^ radicals. The geometries of the optimized 
neutral and charged alanine dipeptides terminated with the —NH^ and —NH2 radicals are 
shown in figuresHand|2lrespectively. In figure^we number different fragmentation channels 
along the polypeptide chain by Roman numerals. These fragmentation channels is a subject 
for investigation in our work. 

The dependence of the total energy of neutral and charged alanine dipeptide terminated 
by the —NH^ radical as a function of distance R between two fragments' for various frag- 
mentation channels is shown in the figures IHHTTl 

Our calculations show that the fragmentation of the alanine dipeptide is a nonmono- 
toneous process. The chain configuration of the neutral dipeptide appears to be not the 
ground state of the molecule as it is clear from the rearangement process occuring in the 
channels I and V. The dipeptide fragmentation process for these channels is exothermic. 
The single ionization of the chain configuration of the dipeptide increases its stability and 
the fragmentation process in this case becomes endothermic. 

The doubly charged alanine dipeptide however becomes unstable because of the Coulomb 
repulsive forces. This fact is illustrated by the energy dependence on R for channel II. 
From this dependence it is seen that the total energy of the system in the final state of 
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FIG. 1: Optimized geometries of neutral (part a), singly charged (part b) and doubly charged (part 
c) alanine dipeptides terminated with the —NH3 radical. 

the fragmentation process decreases as compared to its initial state, which means that the 
doubly charged alanine dipeptide ion terminated with the —NH^ radical is unstable towards 
fission, which should take place if the temperature of the molecule is sufficiently high. 

Note that for many fragmentation channels the energy dependence exhibits a barrier 
calculation of which takes account of the molecular rearrangement. Such rearrangement 
may occur before the actual separation of the daughter fragments begins. In many cases the 
rearrangement leads to the formation of a quasimolecular state (see figures I31 ITT|) . 

Fragmentation of the alanine dipeptide terminated with the —NH2 radical occurs differ- 
ently. Thus the fragmentation of the doubly charged ions of this kind takes place without 
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FIG. 2: Optimized geometries of neutral (part a) and singly charged (part b) alanine dipeptides 
terminated with the —NH2 radical. 

any fission barrier. In the figures IT^lfT!^ we present the dependence of the dipeptide energy 
on R for the fission channel II and V in which the energy release is maximum. 
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FIG. 3: Total energy of neutral alanine dipeptide terminated by the —NH^ radical as a function 
of distance between two fragments for fragmentation channel I (see figure^ calculated with the 
B3LYP method. The energy scale is given in eV, counted from the ground state energy of the 
optimized chain configuration of the neutral alanine dipeptide, Eq "^ = —571.590726 a.u. The 
numbered geometries correspond to the marked points on the energy curve. 
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FIG. 4: Total energy of neutral alanine dipeptide terminated by the —NH^ radical as a function 
of distance between two fragments for fragmentation channel II (see figure^ calculated with the 
B3LYP method. The energy scale is given in eV, counted from the ground state energy of the 
optimized chain configuration of the neutral alanine dipeptide, Eq "^ = —571.590726 a.u. The 
numbered geometries correspond to the marked points on the energy curve. 
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FIG. 5: Total energy of neutral alanine dipeptide terminated by the —NH^ radical as a function 
of distance between two fragments for fragmentation channel III (see figure^ calculated with the 
B3LYP method. The energy scale is given in eV, counted from the ground state energy of the 
optimized chain configuration of the neutral alanine dipeptide, Eq "^ = —571.590726 a.u. The 
numbered geometries correspond to the marked points on the energy curve. 



> 



LU 

I 

LU 



Ala-Ala*° (NH -terminated) [B3LYP/6-31G] 




6 8 10 12 
R (Angstroms) 



18 



1^ ^ 




^♦^•t^ *^ 






X. 



FIG. 6: Total energy of neutral alanine dipeptide terminated by the —NH^ radical as a function 
of distance between two fragments for fragmentation channel V (see figure^ calculated with the 
B3LYP method. The energy scale is given in eV, counted from the ground state energy of the 
optimized chain configuration of the neutral alanine dipeptide, Eq "^ = —571.590726 a.u. The 
numbered geometries correspond to the marked points on the energy curve. 
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FIG. 7: Total energy of singly charged alanine dipeptide terminated by the —NH^^ radical as a 
function of distance between two fragments for fragmentation channel I (see figure^ calculated 
with the B3LYP method. The energy scale is given in eV, counted from the ground state energy 
of the optimized chain configuration of the singly charged alanine dipeptide, E_^-^ ^ = —571.446367 
a.u. The numbered geometries correspond to the marked points on the energy curve. 
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FIG. 8: Total energy of singly charged alanine dipeptide terminated by the —NH^ radical as a 
function of distance between two fragments for fragmentation channel IV (see figure ^) calculated 
with the B3LYP method. The energy scale is given in eV, counted from the ground state energy 
of the optimized chain configuration of the singly charged alanine dipeptide, £'_,_^ ^ = —571.446367 
a.u. The numbered geometries correspond to the marked points on the energy curve. 
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FIG. 9: Total energy of singly charged alanine dipeptide terminated by the —NH^^ radical as a 
function of distance between two fragments for fragmentation channel V (see figure ^ calculated 
with the B3LYP method. The energy scale is given in eV, counted from the ground state energy 
of the optimized chain configuration of the singly charged alanine dipeptide, E_^_-y ^ = —571.446367 
a.u. The numbered geometries correspond to the marked points on the energy curve. 
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FIG. 10: Total energy of doubly charged alanine dipeptide terminated by the —NH-^ radical as a 
function of distance between two fragments for fragmentation channel I (see figure ^) calculated 
with the B3LYP method. The energy scale is given in eV, counted from the ground state energy of 
the optimized chain configuration of the doubly charged alanine dipeptide, E^2 ^ = —570.979569 
a.u. The numbered geometries correspond to the marked points on the energy curve. 
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FIG. 11: Total energy of doubly charged alanine dipeptide terminated by the —NH^ radical as a 
function of distance between two fragments for fragmentation channel II (see figure ^ calculated 
with the B3LYP method. The energy scale is given in eV, counted from the ground state energy of 
the optimized chain configuration of the doubly charged alanine dipeptide, E^2 ^ = —570.979569 
a.u. The numbered geometries correspond to the marked points on the energy curve. 
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FIG. 12: Total energy of doubly charged alanine dipeptide terminated by the —NH2 radical as a 
function of distance between two fragments for fragmentation channel II (see figure EJ calculated 
with the B3LYP method. The energy scale is given in eV, counted from the energy E,2 ^ = 
—570.206214 a.u. The numbered geometries correspond to the marked points on the energy curve. 
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FIG. 13: Total energy of doubly charged alanine dipeptide terminated by the —NH2 radical as a 
function of distance between two fragments for fragmentation channel V (see figure EJ calculated 
with the B3LYP method. The energy scale is given in eV, counted from the energy E,2 ^ = 
—570.206214 a.u. The numbered geometries correspond to the marked points on the energy curve. 
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